The Dynamics of Curved Fronts: Beyond Geometry 
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We derive a new set of kinematic equations for front motion in two-dimensional bistable media. The 
equations generalize the geometric approach by complementing the equation for the front curvature 
with an order parameter equation associated with a nonequilibrium Ising-Bloch bifurcation. The 
resulting equations capture the core structure of spiral waves and spontaneous spiral- wave nucleation. 
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Traveling wave phenomena in reaction-diffusion sys- 
tems often involve sharp interfaces or fronts separat- 
ing different reaction states. The dynamics of two- 
dimensional sharp fronts has been studied successfully us- 
ing a geometric approach jl]-f| . Given a relation between 
the normal velocity of the front and its curvature, the 
geometric theory consists of a closed integro-differential 
equation for the front curvature from which the front line 
shape in the physical plane can be extracted, fnherent 
in this approach is the assumption that the inner front 
structure does not change significantly in time. This as- 
sumption rules out major changes like spontaneous nu- 
cleation of spiral waves along the front. Such phenom- 
ena have been observed recently in numerical simulations 
of model equations describing bistable reaction-diffusion 
systems. Very often the nucleation of spiral waves trig- 
gers spot replication and spiral turbulence 

In this Letter we present a new kinematic approach 
for front motion in two-dimensional bistable media that 
captures spontaneous spiral-wave nucleation along the 
front. A key step in this approach is the consideration 
of a parameter range including a nonequilibrium Ising- 
Bloch (NIB) front bifurcation. This parity breaking bi- 
furcation renders a stationary planar front unstable and 
gives rise to a pair of stable counter-propagating fronts. 
The bifurcation has been found in a number of models, 
including the forced complex Ginzburg-Landau |I(J and 
FitzHugh-Nagumo pd Hl^j e quations, and in experiments 
with chemical reactions |14| and liquid crystals |H| . 

Our kinematic approach consists of three equations: 

• A geometric equation for the front curvature, k: 

If --f+B?* 5 - (i) 

• An equation relating the normal front velocity C„, 
the curvature k, and the order parameter, Co, as- 
sociated with the NIB bifurcation: 

C n =C -DK. (2) 



• An equation for the order parameter: 
dC 

— = (a c - a)C - I3CI + 7« + 7o 

In these equations s is the front arclength, and the critical 
parameter value a c designates the NIB bifurcation point. 
Notice that Co coincides with the planar front velocity 
when k = 0. 

The curvature equation (H) together with the eikonal 
equation (^), where Co is considered constant, constitute 
the geometric approach used in earlier studies H,|[. Re- 
laxing the requirement of constant Co by adding Eq. (||) 
allows for spontaneous local reversal of the direction of 
front propagation. The reversals are accompanied by the 
nucleation of spiral- wave pairs. In the rest of this Letter 
we describe the derivation of Eqs. (||) and (||) for a par- 
ticular model and use these equations to demonstrate a 
mechanism of spontaneous spiral-wave nucleation. 

We consider the FitzHugh-Nagumo model with a dif- 
fusing inhibitor, 

— = e _1 (u -u 3 - v) + S _1 V u , 
dv 

— = u - aw - a + V 2 v , (4) 

where u and v, the activator and the inhibitor, are real 
scalar fields and V 2 is the Laplacian operator in two di- 
mensions. The parameter a\ is chosen so that (|]) de- 
scribes a bistable medium having two stable uniform 
states: an "up" state (u + ,v + ) and a "down" state 
(u— ,V—). Ising and Bloch front solutions connect the two 
uniform states (u± ,v±) as the spatial coordinate normal 
to the front goes from — oo to +oo. The parameter space 
of interest is spanned by e, <5 and do, or alternatively by 
f] = ^/ ^ eS, /i = e/S, and a$. Note the parity symmetry 
(it, v) — > (—it, —v) of (|J) for a = 0. 

The NIB bifurcation line for ao = is shown in Fig. 51 
For (j, <C 1 it is given by <5 = 5p(e) = ?y 2 /e, or r\ = rj c , 
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FIG. 1. The NIB front bifurcation and transverse instabil- 
ity boundaries. The thick line is the NIB bifurcation, (5(e), 
and the dashed lines are the boundaries for transverse insta- 
bility of Ising, Si(e), and Bloch, 5(e), fronts. The thin lines are 
the linear approximations to the transverse instability bound- 
aries near the codimension 3 point, P3. Parameters: a± = 4.0, 
ao = 0. 



where ?y c = 



2\/2g £ 



and q 2 = a x + 1/2 113] . The single sta- 



tionary Ising front that exists for rj > rj c loses stability 
to a pair of counter- propagating Bloch fronts at r\ = rj c . 
Beyond the bifurcation (rj < rj c ) a Bloch front pertain- 
ing to an up state invading a down state coexists with 
another Bloch front pertaining to a down state invading 
an up state. Also shown in Fig. [l] are the transverse in- 
stability boundaries (for ao = 0), 6 = §i(e) = t/rj 2 and 
5 = ^s(e) = rj c /\fe, for Ising and Bloch fronts respec- 
tively. Above these lines, 5 > 5i^b, planar fronts are 
unstable to transverse perturbations |7j,|| . All three lines 
meet at a codimension 3 point P3: e = rj 2 , S = 1, ao = 0. 

The following assumptions are made to derive Eqs. (^) 
and (H): e and S are in the proximity of the codimension 3 
point, P3, with e << 1; the radius of curvature is much 
larger than the front width, that is, k <C q- First, we 
transform to an orthogonal coordinate system (r, s) that 
moves with the front, where r is a coordinate normal to 
the front. Let X(s, t) = (X, Y) denote the position vector 
of the front. More precisely we identify X(s, t) with the 
u = contour line. The relation between the laboratory 
frame (x, y) and the moving frame is 

dY dX 
x = X(s,t) + r— y = Y(s,t)-r—. (5) 
os as 

In the moving frame Eqs. (^) become 

eDu = u — v? — v + \xC,u , 
T>v = u — a\v — ao + Lv , (6) 



where 
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and the curvature is k — ^-%r%r 

os as z 



15Il)sZ~ 101' ^ e cnan g e °f the arclength in time is due 
to stretching and is given by 



ds 

dt 



(7) 



Recalling that jj, = e/S -C 1, we use singular pertur- 
bation theory and distinguish between an inner region 
where ~ 0(1), and outer regions where /if^r ~ 

0(/i). The inner region pertains to the front core where 
the profile of u in the normal direction is steep. Intro- 
ducing a stretched coordinate z = r/^fii and expanding 



u = uq + eu\ + e 2 U2 + ... and v = vq + ev\ + e 2 V2 



we obtain at order unity uo — — tanhz/\/2, 
order e the solvability condition is 
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where Vf — v(Q,s,t) + 0(e 2 ) is the approximately con- 
stant value of the inhibitor v in the narrow [O («//!)] front 
core region. The first term on the right-hand-side of (|§|) 
is identified with the order parameter for the NIB bi- 
furcation: Co = ^^f 11 /- Since the normal velocity is 

C n = — §|) Eq. (||) yields the eikonal equation (||) with 
D = 6- 1 . 

In the outer regions ~ 0(1) and the leading order 
equation for u is u — u — v = 0. The relevant solutions 
are u — u+{v) ~ 1 — v/2 for r < and u — u—(v) ~ 
— 1 — v/2 for r > (assuming a\ is sufficiently large) |12]| . 
To leading order in e we obtain for v the free boundary 
problem 
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and the square brackets denote jumps of the quantities 
inside the brackets across the front at r = 0. 

To solve this free boundary problem we consider a pa- 
rameter range in the immediate vicinity of the P3 point 
in Fig. |l|. In that range the transverse instabilities of the 
fronts involve only small wavenumbers and therefore we 
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can assume weak dependence of v and k on the arclength 
s. In addition, the front speed is small and vanishes at 
P3. This suggests using the speed of a planar Bloch front 
solution, c oc ^Jr\ c — Tj, as a small parameter. The weak 
dependence of v and k on s is achieved by introducing 
the slow length scale S — cs and assuming X = X(5, t). 
This assumption dictates K — c 3 kq where kq ~ 0(1). 
We also introduce a slow time scale T = c 2 t to describe 
deviations from steady front motion. 

Following Ref. p7[ , we solve the free boundary prob- 
lem (^|) by expanding propagating curved front solutions 
as power series in c around the stationary planar Ising 
front 



v(r, S, t, T) = «W (r) + ^ cV n) (r, 5, t, T) , (12) 

n=l 



where t)(°'(r) = (1 - e 9r )/9 2 for r < and «(°)(r) = 
[e~ qr —l)/q for r > 0. Expanding ?7 = r] c — c 2 rji+c 4: ri2+--. 
and using these expansions in (||) produces the set of 
equations 



dt ' * <9r 2 
where and p( 2 ) are as in Ref. Jl 



= -p {n \ n=l,2,3,. 



(13) 



and p 
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In ( jl4| ) we assumed a = c a 0Q where a Q0 ~ O(l), and 
recall that k = k/c 3 . Notice that ^ ~ C(r)i and 
P2 contributes only at orders higher than c 3 . We solve 
Eq. (|l3|) using the asymptotic behavior of an appropriate 
Green's function as in Ref. |l7]] . The results for n = 1 , 2 
remain unchanged and give the front bifurcation point 
rjc = - JL 3 . The solution of jl^ ) with n = 3 yields the 
compatibility condition 
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or expressing the slow time and arclength derivatives in 
terms of the fast variables t, s and using (0) and ([l2"|), 
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Equation (|T^) coincides with (||) once we make the follow- 
ing identifications: C = ~^Vf, « = ^1, «c = 

/3 = 1/6, 7 = a c (l - S^ 1 ), and 70 = 2a c qa a . 

Equation (||) reproduces the NIB bifurcation for planar 
fronts: setting k — and ao = we find the Ising front 
branch Co = and the two Bloch front branches Co = 
±y/(a c — a)/ [3. To test whether Eqs. (0)-(U) also capture 
the transverse instabilities we check the linear stability 
of planar front solutions near the P3 point in the ao = 
plane. Let Co = Cq + Co exp(ert + iQs) and n = k° + 
Rexp(at + iQs) where (Cq, k°) = (0, 0) for the Ising front 
and (C$,k°) = (±v/(a c - a)//3,0) for the Bloch fronts. 
Inserting these forms in (^) gives the following transverse 
instability lines, linearized around 5=1: 

Ising : e = r? 2 <5 , Bloch : e = rj 2 (3 - 25) . 

These lines are displayed in Fig. [I] (thin lines). To linear 
order around the P3 point they coincide with the exact 
transverse instability lines. 

As a first application of the kinematic equations (0)-((|) 
consider a "front" solution connecting the planar Bloch 
front, Co = Cq , k = 0, at s = —00 with the pla- 



nar Bloch front, Co = C 



0, at 



-00, where 



C, 



ct)//3, and we have assumed a symmet- 
ric model, ao = or 70 = 0. Fig. ||a shows such a solution 
obtained by numerically integrating (0)-(§). As demon- 
strated in Fig. 2]b this front solution of the kinematic 
equations ((l])-(^) represents a spiral-wave solution of the 
FitzHugh-Nagumo model (^) . Unlike the geometrical ap- 
proach [^],Q the spiral core is naturally captured by the 
new kinematic equations. 



0.1 



0.0 



-0.1 



C A 


(a) 








200 



600 



-100 100 
s X 

FIG. 2. A front solution to the kinematic equations 
(a) The order parameter Co and the curvature k along the ar- 
clength s. (b) In the x — y plane the front solution corresponds 
to a rotating spiral wave. The shaded (light) region corre- 
sponds to an up (down) state. Parameters: ai — 4, ao = 0, 
e = 0.01234, 5 = 1.0. 

A second application of the kinematic equations is 
the study of spontaneous spiral-wave nucleation. Spiral- 
wave nucleation, induced by a transverse instability, has 
been previously observed in direct simulations of (m M. 
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Figs. |]a-d show the time evolution of a solution to the 
Co — k equations representing a planar front near the 
NIB bifurcation and beyond the transverse instability 
boundary. The initial front pertains to an up state in- 
vading a down state (Co > 0). The transverse instabil- 
ity causes a small dent on the front to grow (Fig. ||b). 
The negative curvature then triggers the nucleation of a 
region along the arclength where the propagation direc- 
tion is reversed (Co < 0) (Fig. ||c). The pair of fronts 
in the kinematic equations that bound this region cor- 
respond to a pair of counter-rotating spiral waves in the 
FitzHugh-Nagumo equations (Fig. |3|d). With this ap- 
proach, the two-dimensional spiral-wave nucleation prob- 
lem is reduced to the considerably simpler problem of 
domain, or droplet, nucleation in one dimension |18|. 
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FIG. 3. Nucleation of a spiral-wave pair in the kinematic 
equations (0)-(^) Left column: the Co(s) and «(s) profiles. 
Right column: the front line shape in the x — y plane. Pa- 
rameters: ai = 4, a = -0.0001, e = 0.0115, S = 1.063. 
(a)-(d) are at t = 0, 116, 136, 142. 

We have derived kinematic equations for front motion 
in two-dimensional bistable systems near a NIB bifur- 
cation. The equations generalize earlier derivations and 
capture both the core structure of spiral waves and the 
dynamic process of spiral-wave nucleation. Further in- 
vestigation is needed to determine if other features of 
spiral waves, like the meander instability JhJ, are cap- 
tured by the equations. Note that front interaction ef- 



fects are excluded by the choice of the boundary condi- 
tions, i;(=foo, s, £) = v±, in Eqs. (||). Such interactions 
are not significant for the initial stages of spiral wave 
nucleation or for the symmetric (or nearly symmetric) 
low curvature spirals studied in this Letter. Front in- 
teractions, however, do become significant when highly 
curved spirals develop, and might play an essential role 
in the meander instability. Goldstein et al. |^0| have re- 
cently studied front interactions in the fast inhibitor limit 
(e 1) where stationary patterns prevail. A combina- 
tion of the approaches used in these two complementary 
studies may prove useful in establishing a theory of spiral 
waves of wider validity. 
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